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Abstract: The Low Frequency Instrument (LFI) is an array of cryogenically cooled radiometers on 
board the Planck satellite, designed to measure the temperature and polarization anisotropics of the 
cosmic microwave backgrond (CMB) at 30, 44 and 70 GHz. The thermal requirements of the LFI, 
and in particular the stringent limits to acceptable thermal fluctuations in the 20 K focal plane, are a 
critical element to achieve the instrument scientific performance. Thermal tests were carried out as 
part of the on-ground calibration campaign at various stages of instrument integration. In this paper 
we describe the results and analysis of the tests on the LFI flight model (FM) performed at Thales 
Laboratories in Milan (Italy) during 2006, with the purpose of experimentally sampling the thermal 
transfer functions and consequently validating the numerical thermal model describing the dynamic 
response of the LFI focal plane. This model has been used extensively to assess the ability of LFI 
to achieve its scientific goals: its validation is therefore extremely important in the context of the 
Planck mission. Our analysis shows that the measured thermal properties of the instrument show 
a thermal damping level better than predicted, therefore further reducing the expected systematic 
effect induced in the LFI maps. We then propose an explanation of the increased damping in terms 
of non-ideal thermal contacts. 
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1. Introduction 

Planck is a space mission of the European Space Agency (ESA) whose main objective is to image 
the temperature and polarization anisotropics of the Cosmic Microwave Background (CMB) with 
unprecedented sensitivity, angular resolution, sky coverage and frequency coverage [14]. Two 
instruments share the focal plane of the Planck 1.5 m telescope: HFI (High Frequency Instrument, 
[6]), an array of 52 bolometers; and LFI (Low Frequency Instrument, [2]), an array of 22 pseudo- 
correlation radiometers. Both HFI and LFI need to be cooled to cryogenic temperatures (0. 1 mK for 
the HFI bolometers and 20 K for the LFI radiometers) in order to meet their scientific requirements. 
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Figure 1. Left: schematics of the Planck satellite. The warm service module (~ 300 K) is thermally 
decoupled from the focal plane and the telescope by means of three conical radiators. The LFT focal plane is 
cooled to 20 K by a hydrogen sorption cooler which also acts as a pre-cooling stage for the HFT 4 K cooler 
that cools the HFI focal plane and the LFI reference load system [17]. Right: detailed view of the LFI 
structure. On top, the cold Focal Plane Unit (FPU) with both the LFI and HFI feed horns is shown. A set of 
waveguides cormect the FPU with the warm (300 K) Back End Unit (BEU), shown at the bottom. 

Planck is the first space mission devoted to the measurement of CMB anisotropics that uses 
cryogenic instruments, and it presents extreme challenges in its thermal design. The thermal system 
uses a mixed passive/active approach to reach the required working temperatures. Three thermal 
shields radiate heat and decouple the warm service module (~ 300 K) from the telescope, passively 
cooled to <50K, and the focal instruments. Active cooling is provided by three cryocoolers: a 
hydrogen sorption cooler (18-20 K), a Stirling cooler (4 K) and a ^He-'^He diluition cooler (0.1 K). 

Within such a complex cryogenic system, thermal effects are expected to be the most important 
source of systematics. Temperature instabilities in the detectors can degrade the scientific output 
of Planck, which must be able to measure fluctuations in the CMB signal with a sensitivity of 
5T IT ~ 10"^ per pixel on the final maps. For this reason, one of the main driving requirements 
in the design of Planck has been to maximize the damping of thermal fluctuations propagating 
through the two instruments, and an extensive set of thermal tests has been performed at both 
instrument and satellite level in order to verify that the thermal stabiUty requirements are satisfied. 

The LFI focal plane, cooled to 20 K, includes front end passive components (corrugated feed 
horns and orthomode transducers), hybrid couplers, and state-of-the-art HEMT (High Electron 
Mobility Transistor) low noise cryogenic ampUfiers providing ~ 30 dB of amplification. The ther- 
mal susceptibility of the front end components, in particular the thermal coupling to the gain of 
the HEMT amplifiers and to the insertion loss of the front-end passive components, impose very 
stringent temperature stability requirements (~ yuK level) on the 20 K stage. 
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One of the most important systematic effects in the Planck measurements is the propagation 
of fluctuations from the 18-20K sorption cooler [7], whose stability at the cold end is limited by 
the thermal cycling of the sorbent compressors [3, 12]. In fact, fluctuations in the LFI 20 K stage 
introduce a potentially serious systematic effect in the LFI science as they may mimic brightness 
changes in the beam as the satellite scans through the sky [15]. While the pseudo-correlation 
design of the LFI receivers suppress to first order thermal fluctuations [9], residual effects must 
be damped to extremely low levels (see table p. The combination of radiometer susceptibility, 
expected fluctuations in the 20 K stage, and thermal damping from the instrument were the key 
factors in the instrument design to ensure adequate stability. 

To quantify thermal damping factors in the focal plane, a numerical thermal model was devel- 
oped by Thales Alenia Space. The model has been used extensively since the design phases (1) to 
study how thermal systematic effects propagate to the radiometer front ends, and from the results of 
this analysis (2) to estimate the stability requirements needed for Planck/LFI to meet its scientific 
objectives. 

In this paper we report the measurements of the propagation of fluctuations in the LFI flight 
model focal plane after the LFI was assembled in the Thales Alenia Space laboratories in Milan 
(Italy). The first product of our analysis is the sampling of the dynamical thermal transfer function 
between the focal plane cold end and a set of points on the focal plane itself at three frequencies. 
We then use such measurements to validate the numerical thermal model of the focal plane, thus 
confirming (1) the validity of the temperature stability requirements on the Planck/LFI focal plane, 
which have driven the design of Planck, and (2) the validity of a number of works that have used the 
estimates of the model [9, 16, 10, 1]. Finally, we use the measured transfer functions to estimate the 
improvement in the stability of the receivers over the LFI requirements and provide an explanation 
of such improvements in terms of non-ideal thermal contacts. 

The outline of this article is as follows: Sec. describes the Planck/LFI instrument, the cooling 
system of Planck and the thermal model used for the characterization of the LFI focal plane. 
Sec. describes the extraction of the focal plane thermal transfer functions and the validation 
of the thermal model of the focal plane. Sec. ^ shows how much the measured transfer functions 
allow to relax the reqmirement on the thermal stability of the radiometers. In Sec. ^ we propose our 
explanation for the discrepancies between the model and the measurements (i.e. better damping) 
in terms of contact resistance between the cold end of the focal plane and the focal plane itself. 
Finally, in Sec. ^ we report the conclusions of our work. 

2. The Planck/LFI Instrument 
2.1 Overview of Planck/LFI 

The Planck/LFI instrument [2] is an aiTay of 22 pseudo-correlation radiometric receivers centered 
at 30, 44 and 70 GHz and cooled to ~ 20 K by a vibrationless hydrogen sorption cooler [12]. Each 
LFI receiver measures the sky signal (~ 2.7 K) by comparing it with the signal coming from a stable 
reference load (~ 4.5 K). This load is thermally connected to the 4 K HFI shield. 

In order to minimise power dissipation in the focal plane, the LFI Radiometric Array Assembly 
(RAA) is split into two subassemblies (see fig. p: the Front End Unit (FEU), mounted on the focal 
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Table 1. Total error budget allocated to thermal fluctuations in the LFl instrument [2]. All the numbers 
are in juK. The first three columns contain the maximum error for high-frequency fluctuations (h.f.), i.e. 
> 1/60 Hz. The last two columns contain the error budget for low-frequency periodic fluctuations and spin- 
synchronous variations. The BEU (Back-End Unit) comprises the warm part of the radiometers (~300 K), the 
DAE (Digital Acquisition Electronics) is the back-end module devoted to the digitization of the radiometric 
signal. 



Source 


H.f. @ 30 GHz 


H.f.@44GHz 


H.f. @ 70 GHz 


Periodic [pK] 


Spin-s. [juK] 


Focal plane 


14.8 


20.5 


36.4 


0.9 


0.45 


Waveguides 


14.8 


20.5 


36.4 


0.4 


0.4 


BEU 


14.8 


20.5 


36.4 


0.4 


0.4 


DAE 


14.8 


20.5 


36.4 


0.4 


0.4 


Total 


29.6 


41.0 


72.8 


1.1 


0.8 



plane and actively cooled at 20 K, and the Back End Unit (BEU) on top of the Planck Service 
Module at about 300 K. The FEU is connected to the BEU by means of waveguides (WGs), which 
carry the microwave signals and by means of the cryo harness that carries the bias currents for the 
front-end active components. Both the WGs and the cryo harness are connected to three passive 
radiators (V-grooves) that dissipate the heat coming from the BEU. 

The goal of carefully controlling thermal systematic effects has driven the design of LFI, which 
has been optimized in order to maximize the stability to thermal fluctuations, especially those 
propagating through the focal plane (see fig. 0). In order to improve the focal plane stability, a 
Temperature Stabilization Assembly (TSA) has been implemented in order to reduce temperature 
fluctuations at the sorption cooler cold end by one order of magnitude. Such design and optimiza- 
tions have been made by means of stability estimates (described in Sec. ^ based on a numerical 
thermal model of the focal plane. The latter is the main topic of this article. 

In order to verify that the thermal performance of LFI is compliant with the requirements, 
the instrument has undergone an extensive thermal analysis and test campaign to evaluate if the 
measured thermal damping is comparable (or better) to thermal model estimates. This is of capital 
importance: temperature fluctuations propagating to the amphfiers can induce gain and insertion 
loss changes that in turns produce systematic variations in the output antenna temperature [15]. 
Therefore, any uncontrolled fluctuation at this level might compromise the whole scientific output 
of the instrument. 

2.2 The Thermal Model of the LFI Focal Plane 

A detailed thermal model of the focal plane was developed during the design of LFI by Thales/Alenia 
Space using the ESATAN/ESARAD [4] numerical thermal analysis software. This model has two 
purposes: (1) to characterize the thermal steady state of each part of the instrument, and (2) to 
quantify dynamical thermal transfer functions between two physical points of the instrument. Both 
derive from the need to allow Planck to control thermal systematic effects at the sub-/iK level (see 
table [I]). In the specific case of the focal plane stability requirements [2], these call for periodic 
fluctuations to be less than ±0.9 /uK per pixel in the final maps and for spin-synchronous (s.-s.) 
variations to be below ±0.45 /uK per pixel. The tighter requirement on s.-s. fluctuations is due 
to the fact that their period is equal to the spin period of the satellite, i.e. 60 s, and therefore the 
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Figure 2. The LFI focal plane during the integration of the flight instrument in the Thales Alenia Space 
laboratories (Milan, Italy). Note the 1 1 feed horns (2 centered at 30 GHz, 3 at 44 GHz and 6 at 70 GHz). 
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Figure 3. Placement of the temperature sensors in the LFI focal plane. One of the sensors (TS5R) has been 
placed on the flange that holds feed horn #28 (left). Two sorption coolers are used in Planck; the position of 
the cold end used during the LFI RAA test campaign (LVHX2) is shown with a black point on the bottom 
left. All the sensors were included in our analysis, with the exception of TS2L and TS5L (see text). 



associated errors in the maps are not reduced by the redundancy of scanning the same sky circle 60 
times. This imposes a high degree of temperature stability at the level of the 20 K sorption cooler 
cold end (< 100 mK peak-to-peak). 

Due to the criticality in measuring the focal plane thermal stability, 12 silicon diode thermome- 
ters are placed on the LFI focal plane (see fig. 0). All sensors have dedicated calibration curves 
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and each exhibits an accuracy of ~ 20 mK. Five of them have been calibrated in the range 14-26.5 
K and the associated readout electronics allow a sensitivity of 0.9 mK. The others have been cal- 
ibrated over the 15-90 K range, with two main sensitivity ranges: 13.3 mK from 25 to 90 K and 
1.4 mK from 15 to 25 K. The intrinsic noise of each thermometer is always below the calibration 
uncertainties. The characteristics of the thermometers are not sufficient to detect the very slight 
temperature changes expected during flight, which are estimated by the thermal model to be of the 
order of 1 mK peak-to-peak. However, they have been used extensively in this work to measure 
the response to induced temperature fluctuations on the FPU (much greater than 1 mK) and will be 
used during flight to monitor any unexpected change in the temperature of the focal plane. 

3. Estimation of the Focal Plane Transfer Functions and Dynamic Validation of the 
Numerical Model 

After the integration of the LFT flight instrument, we verified the ESATAN thermal model of the 
focal plane by forcing thermal fluctuations at the sorption cooler cold end in cryogenic conditions, 
measuring the induced fluctuations on the thermometers placed on the focal plane (see fig. ||) and 
then comparing the results with the estimates of the thermal model. Such analysis has allowed us 
to verify the compliancy of the LFI focal plane temperature stability with scientific requirements, 
as well as to experimentally sample the thermal transfer function for a set of frequencies. 

3.1 Test Methodology 

To perform the measurements, we induced a sinusoidal temperature change of fixed frequency v, 
amplitude ATq and phase (po (the input fluctuation) at the cold end and measured the response at 
each of the focal plane thermometers. 

In principle, any temperature profile can be induced and thermal transfer functions can be 
derived through a Fourier transform, which allows multiple frequencies to be studied at the same 
time. However, preliminary studies during the calibration of the LFI Qualification Model (QM) 
in 2005 proved that variations with a strongly peaked spectrum are easier to analyze and produce 
more accurate results. 

From the time stream of temperatures (one for each thermometer) we calculate the amplitude 
AT and phase (p of the sinusoids and compared them with the input fluctuation. We estimate the 
y factor as the ratio of the two amplitudes AT/ ATq and the phase shift Aip as the difference of 
the two phases (p-(po- By repeating this test with different frequencies v for the input fluctuation, 
we reconstruct the profile of y(v) and Aip{v) as a function of the frequency. The complex quantity 
y(v)exp(/^(v)) is called the thermal transfer function between the cold end and the thermometer' 

During the LFI FM cryogenic tests in the Thales Alenia Space laboratories we performed three 
separate injections of sinusoidal fluctuations with frequencies O.lSmHz, 0.55 mHz and 1.40 mHz 
(corresponding to periods of 12, 30 and 60 minutes) in the cold end. This was done by inducing a 
time-varying power of sinusoidal shape through a resistor fixed near the cold end. 

To derive the transfer function from the thermal model, we apply a sinusoidal temperature 
profile as a boundary condition and calculated the predicted damping and phase shift. 

'it is easy to show that if ye"'' is the transfer function from point 1 to point 2, then the t.f. from point 2 to point 1 is 
y-lg-up gjj^jj jj^jg relation in Sec. pi 
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Figure 4. Example of a fluctuation induced at the cold end of the focal plane and the measured temperature 
response at one of the thermometers (TS5R, placed on the flange of feed horn 28, see Figure 0) during 
the LFI RAA tests. Above: Temperature profiles at the cold end and at the flange. The amplitude of the 
fluctuation at the flange is smaller than the one at the cold end because of the thermal path. Below: Spectral 
profile of the two temperature profiles. Each profile has a peak centered at v = 0.55 mHz, corresponding 
to a period of 30 minutes. (The 1 // component is due to the cold end and not to the intrinsic noise of the 
thermometers.) 



3.2 Numerical Methods used in the Analysis 

The standard method used to estimate transfer functions is based on a straightforward application of 
the Fourier Transform. This kind of analysis has the advantage of being quite simple to implement 
and produces solid results but, on the other hand, it does not provide an easy way to quantify 
errors in the estimation of the transfer functions. Such eiTors are generated mainly by the slow 
temperature drifts, the signal quantization induced by the digital acquisition board (see fig. ^ and 
the intrinsic noise of the thermometers. For this reason, we have developed three different analytical 
methods that can both reduce the impact of such effects on the final result and quantify the errors 
of the estimation: (1) a Fourier method which uses a jackknife-hke test to derive the error on the 
estimates of y and A^, (2) a method which directly extracts y and Aip from time-domain data and 
(3) a method using non-linear fitting algorithms. The three methods are explained in Appendix 0. 

To derive our best estimate for y and Aip we simply picked the estimate with the lowest relative 
error (e.g. 6y/y) among the one produced by the three methods. In general, the fitting method 
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Figure 5. Temperature measured by the thermometer near the LFI feed horn #25. The temperature sinusoid 
apphed at the cold end had a frequency of 0.55 mHz. Note the quantization induced by the digital ther- 
mometer and the slow thermal drift (longer than one period of the sinusoid). Any method used to extract the 
amplitude and phase of the sinusoid must be able to limit the impact of both effects on the calculations. 

has produced the best results for y while the Fourier method has outperformed the others in the 
determination of A^o. 

Calculating y and from the numerical thermal model is straightforward, because these tem- 
perature streams suffer neither quantization errors nor drifts. Therefore, for the sake of simplicity 
we have analyzed numerical data only using the time-domain method. 

3.3 Discussion of tlie Results 

The purpose of the validation tests on the LFI focal plane model is to verify that the measured 
damping level of thermal fluctuations did not overcome the estimate of the numerical model, since 
the latter has been verified to be compatible with the scientific requirements of LFI. 

Of the 12 thermometers in the focal plane, we have chosen to exclude from our analysis TS2L 
and TS5L because of a few problems occurred with them during the acquisition. 

The results for the analysis of the 10 sensors for each of the three frequencies are reported in 
tables H, ^ and ^ at the end of the article. The three methods show an excellent agreement (within a 
few percent) for their estimates and errors on y and y. We get larger error bars for A^ when using 
the fitting method (see table however, the three methods produce estimates which are always 
within 2o", therefore confirming the validity of each method. 

The comparison between these values and the estimates of the numerical model is shown in 
figures and [TT]. In these plots y is always smaller in the experimental Une than in the numerical 
estimates, while the contrary appUes to A(/j. Both effects were somewhat expected, because the 
thermal model cannot take into account unknown contact resistance effects. This topic will be 
addressed in Sec. 0. 
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Parameter 


Value 


iNotes 


Mli-OMT 


(J. 25 dB 


Measured at room temperature 




0.25 dB 


Measured at room temperature 


^sky 


3.7 K 


2.7 K from the sky plus 1 K from the telescope 




4.5 K 


Estimated value 


jFEM 
phys 


26 K 


From table 3 in [11] 


Tn 


10.5 K 


From table 1 


in [8] (mean between 28-MO and Ml) 


dG 

at 

phys 

r 


-0.03 dB 


From table 11 in [15] 


0.15 dB 


From table 11 in [15] 


0.95 


Estimated value 


7-FEM 


1.05x10-2 


From eq. 


4.1 





Table 2. List of parameters used in eq. ^[T] to estimate the radiometric transfer function for radiometer 
#28-M (30 GHz). 

4. Impact of Instabilities on LFI Maps 

In this section we use the measured transfer functions to quantify how the better damping (i.e. 
smaller values for y) will reduce the error caused by sorption cooler instabilities in the LFI maps. 
To do this, we need to simulate the propagation of fluctuations at the sorption cooler cold end to 
the production of the maps. The stpng involved in this calculations are the following (taken from 
[10]): 

1 . the temperature fluctuation propagates from the cold end to the radiometer in a way described 
by the thermal transfer functions discussed in this paper; 

2. the fluctuation in the radiometer temperature induces a fluctuation in its output which is 
proportional to the temperature itself [15]; 

3. in the production of the map, multiple passes over the same pixel are averaged, thus further 
reducing the error caused by non-spin synchronous fluctuations; 

4. destriping techniques are going to be applied to the Planck maps in order to reduce the impact 
of long-term drifts, which include sorption cooler fluctuations as well. 

4.1 Assumptions Used in the Analysis 

We have studied the case for radiometer #28-M (30 GHz), as horn #28 is the easiest to use for our 
study because of the temperature sensor mounted on its flange (sensor TS5R, fig. To carry out 
the calculations, we use the following assumptions: 

• The radiometric response rj™ is defined such that the temperature change on the map 6T^ap 
due to a variation of the physical temperature ^Jphys is Jr^ap = x (Jrphys. We use eq. 2.5 
in [15] to estimate analitically, with the simplifying assumptions of a perfect match in 
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the two gains (G :- Gpv = Gf2) and noise temperatures {T„ := T„fi = T^fi)' 



rP™ = Lffi.oMTx||l--^|-r|l- — 1 + 



with 



/ 



1 \ / 1 



-f-a-OMT / \ ^4K / 



r ~ ~ 1 dG dT„ \ 

[Ts^y + Tn - r{T4K + Tn)\ + ( 1 - r) ^ 

phys phys ^ 



Tsky I 1 



i sky = + 1 - Tpi^ys 

i-fh-OMT \ i-fh-OMT ' 



nFEM 



'4K \ MK 



(4.1) 



The Ust of parameters used in the estimation of the radiometric t.f. are listed in table 0. 

The measurement redundancy is estimated using an analytical model instead of considering 
the full details of the Planck scanning strategy. This model only estimates the redundancy of 
the pixels on the ecliptic plane, which is the worst case. The transfer function is the following 
(taken from [10], eq. 3): 



sin(:TAfy/yspin) 



(4.2) 



sin(;Ty/Vspin) 

where A'^ is the number of scan circles over which the average of each pixel is computed (in 
this context we use A'^ = 24 x 60 = 1440, corresponding to a 24-hour scanning period over the 
same circle) and v^pm = 1/60 Hz is the Planck spin frequency. 

The impact of destriping on the maps is estimated via a transfer function, instead of applying 
the full algorithm to the map data. We use the following transfer function (taken from [10], 
fig. 5): 



Tf'''{v) = -+b, (4.3) 

J y 



a 



with a - 0.1067 Hz and b - 1.7992. 



With such assumptions, a fluctuation amplitude 5rcoid-end with frequency y at the cold end 
induces an error in the map equal to 

<5r,,ap - (Tj^'^v) rp^iv) Tf^ r(v)) X 5r,oM.end, (4.4) 

with y(y) being the thermal transfer function between the cold end and the radiometer. As said 
above, in our case we shall use the values of y reported in table ^ for sensor TS5R. 

4.2 Study of Sinusoidal Fluctuations 

In this paragraph we use the algorithms and formulae described above to study how the fluctuation 
amplitudes reported in table [I] change when switching from the numerical to the measured transfer 
functions. Our approach follows these stpng: 
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V 


<5^cold end 


^ map 


ry^meas 
^ map 


Improvement 


0.25 mHz 


0.339 mK 


0.900 /iK 


0.526 juK 


41% 


0.60 mHz 


0.244 mK 


0.900 /iK 


0.549 /iK 


39% 


1.50 mHz 


0.257 mK 


0.900 /iK 


0.539 /iK 


40% 



Table 3. Difference between the peak-to-peak fluctuation in the LFl temperature maps estimated using the 
numerical (6T^™) and the measured (ST^^p'^) thermal t.f. for three frequencies. The 0.25 mHz and 1.5 mHz 
frequencies are respectively the average frequency of each of the six sorption cooler compressor elements 
cycle and of the overall compressor cycle, and they are therefore of considerable importance in the analysis 
of LFl systematics. The fluctuation amplitude at the cold end ^rcoidend has been chosen so that a 24-hour 
observation leads to 6T^^ - 0.900/iK. The calculation has been done for radiometer #28-M (30 GHz). 

1 . By inverting the whole process, we derive a fluctuation amplitude at the cold end that pro- 
duces an error in the map equal to the number reported in table [I] when the numerical t.f. is 
used; 

2. We study the error on the map caused by the same fluctuation when the measured t.f. is 
applied. 

We have chosen not to do this compaiison for spin-synchronous and high-frequency fluctuations, 
as our measurements do not cover such frequency ranges. Instead we concentrated on a range of 
frequencies similar to the ones measured during the instrument tests. 

Table reports the results of our calculations: the reduction in the value of j{v) from the 
numerical to the measured numbers (between 39% and 44%, see fig. [T^ has lead to a comparable 
improvement in the magnitude of the error in the map, between 39% and 41%. 

5. Estimation of the Impact of Contact Resistances in the Measurements 

The discrepancies between the model and the experimental measurements in general can be ex- 
plained by the presence of non-ideal phenomena that are often very hard, if not impossible, to 
simulate. Among these phenomena, in a thermal system, the contact resistance is usually one of 
the most important. In this section we shall expand our idea proposed in Sec. that the differ- 
ences between the model and the data are due to the presence of unknown contact resistances in 
the system and in particular at the main interface between the sorption cooler cold end and the LFl 
focal plane. 

A contact resistance is a constant that characterizes a non-ideal thermal contact between two 
bodies. It is expressed as the ratio between the temperature jump Ar at the interface and the flux 
per unit area q across the interface itself: 

Ar 

Rc = —. (5.1) 

As shown in Appendix 0, when a sinusoidal fluctuation propagates through the interface (the case 
of interest in the context of this paper) there is a sharp decrease in the value of the fluctuation 
amplitude y and an increase in the value of the phase shift A^ (see fig. Note that this is in 
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agreement with our comparison between estimated and measured transfer functions, as the mea- 
surements always showed greater damping (smaller y, see fig. [T^) and larger phase shift (larger A^, 
see fig. \n\ ). 

The experimental measurement of contact resistances is a complex task which requires ded- 
icated instrumental set-ups and is therefore not easily applicable in the context of the Planck/LFI 
focal plane thermal tests, where a limited number of temperature sensors are available and the in- 
volved objects have complex shapes. From eq. pTTj in order to estimate Rc the temperature change 
at the interface and the heat flux flowing through the interface must be estimated. Both quantities 
are not measurable directly but must be extrapolated from a number of temperatures measurements 
near the interface^. Although it is sometimes possible to estimate Rc from transient temperature 
measurements (see e.g. [5]), the most common set-up uses bodies in thermal equilibrium in the 
one-dimensional approximation, since this configuration is easily modeled mathematically (see 
Appendix |B.1D . In any case, for the mathematical problem to be solvable it is required to measure 
the temperature at a number of points near both sides of the interface: a quick glance at fig. 
reveals that this is very hard to accomplish with the Planck/LFI focal plane thermometers. 

To increase our confidence in explaining the discrepancies observed in the measured vs. esti- 
mated transfer functions, we have therefore chosen to perform a statistical analysis of such discrep- 
ancies. We have estimated all the possible transfer functions between two sensors of the focal plane 
or the sensor at the cold end (both numerically and from the same experimental data discussed in 
the previous sections) and have divided them into two groups: 

Focal plane t.f.: These are the transfer functions between two thermometers on the focal plane, 
i.e. pairs of the thermometers shown in fig. ^ Since there are 10 thermometers available, this 
leads to 10 x 9 = 90 transfer functions. For the sake of simplicity we have estimated the t.f. 
using the time-domain method (see Appendix ^ instead of applying the three methods and 
picking the best estimate. 

Cold end t.f.: These are the transfer functions between the cold end sensor and one thermometer 
on the focal plane. They have been derived and discussed in Sec. 0. 

Since non-ideal thermal contacts are not considered by the numerical model, if we expect that an 
important non-ideal contact be at the cold end/focal plane interface then the discrepances between 
the numerical and the measured t.f. should be larger for the cold end transfer functions than for the 
focal plane t.f. We quantify the discrepance between two different estimates of a t.f. (i.e. numerical 
vs. measured) in terms of two quantities: 



where (ymeas(^)> A(^meas(v)) is the measured t.f. and (Ai^num(>'). Ai^num(v)) is the numerical estimate. 
(With this definition of dy and d^, the two quantities do not change when one exchanges the two 

"An additional complication is that the mathematical problem is ill-posed and requires the use of inverse method, 
which often suffer from numerical instabilities. See [13]. 



6y{v) - 2 



T lymeas(v)-ynum(v)l 

rmeas(v)-l-rnum(v) ' 
\Aip 

meas 

(v) - Aip 

num (V)l 



(5.2) 



(5.3) 
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sensors, i.e. 5^^^ = dy~^^ and 5^"^^ = S^~*K) In principle each quantity is able to estimate the 
discrepancy between two t.f. However, we expect the presence of the contact resistance to have a 
greater effect in 6y than in 6^ (see Sec. |B.2| in the appendix and expecially fig. ^: compare e.g. the 
jump in y{v) and in A(p{v) for Rc = 1). 

The result of our analysis is shown in fig. ^ which shows the points with coordinates ((5y,5y,) 
for each of the three frequencies and each of the focal plane t.f. (white circles) and the cold end t.f. 
(black circles). The convex hulls'' of the two sets of points are shown as well. It is evident that if we 
consider the 6y parameter (abscissa), then discrepancies in cold-end t.f. are systematically higher 
than those in focal plane t.f.: this is exactly what we were expecting from our hypotheses, and 
thus confirms the fact that the differences between the measurements and the numerical estimates 
in figures [T^ and [TT] are likely to be due to a non-ideality at the SCS cold end/LFI focal plane 
discontinuity. 

6. Conclusions 

We have presented the results of our analysis of the thermal performance of the focal plane of the 
Planck/LFI instrument. The purpose of this work is twofold: (1) to experimentally measure the 
thermal transfer function between the 20 K focal plane cold end and a set of points in the focal 
plane itself, and (2) to use such measurements to validate the numerical thermal model of the focal 
plane. The latter point is of capital importance for LFI, since the design of the instrument has been 
driven by the thermal stability estimates produced by the model itself. 

We have estimated the focal plane thermal transfer functions by inducing a sinusoidal tem- 
perature variation near the cold end while measuring the impact on the thermal stability on other 
points of the focal plane. We have applied three different analysis methods which have provided 
results in good agreement, especially in estimating the amplitude reduction (the y factor). 

Our analysis has compared measured transfer functions with simulations produced via the 
numerical thermal model of the LFI focal plane. We have found that in every case the thermal 
mass of the instrument is able to damp fluctuations better than predicted by the model, probably 
because of a non-ideal contact between the focal plane and the cold end. Therefore, this analysis 
not only confirms the LFI thermal performance compliancy to the expected sensitivity per pixel 
level of AT/T ~ 10"^ in the final maps, but in some cases estimates an improvement of ~40% in 
the peak-to-peak effect of long-term focal plane temperature fluctuations on the maps produced by 
LFI. 

The next step of our study is to assess the impact of the residual effect of such fluctuations 
on the science, as well as to characterize the impact of fluctuations on the LFI focal plane induced 
by thermal instabilities originating within HFI, the other instrument on board of Planck. Both 
activities will be possible by exploiting data taken during integrated satellite tests at cryogenic 
conditions on ground and during the early flight phases, currently being analyzed. 
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Figure 6. Bivariate scatter plot of the quantities 5y and d^p (equationsp.2|and|5.3p for 10 out of 12 focal plane 
thermometers. They have been grouped according to the three frequency measured during the instrument 
tests. White points have been calculated using the transfer functions between sensors on the focal plane 
(10x9/2 — 45 points), while black points have been calculated using the t.f. between the cold end sensor 
and each of the 10 focal plane thermometers used in this work. The white and grey polygons are the convex 
hulls of the two point sets. The discussion of these plots is provided in the text. 
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A. Analysis IVIethods 

Here we provide a short description of the three analysis methods used to extract the cold-end 
thermal transfer functions to the focal plane thermometers. 

A.l Fouriermethod 

The Discrete Fourier Transform (DFT) of the temperature streams at the two points is calculated 
after the two datasets have been interpolated to the same sampling frequency. The value of y is 
given by the ratio in the height of the peaks in the two transforms, while the estimate of the phase 
delay A(p is given by the difference in the phases of the two peaks. 

As said above, getting a meaningful estimate of the error of y and A(p from this method is not 
trivial. We calculated v, y and Acp over A'^ partially overlapping time windows, thus obtaining N 
estimates for each quantity. We then considered our best estimate and its error as the mean value 
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Figure 7. Layout of the time windows used in the Fourier method. Within a time window of length T, we 
consider intervals of length (N-l)/NxT. Each interval j -I...N starts eit (j - 1)/{N - l)xT /N. 

and the standard deviation of the set, respectively^. A sketch showing how the A'^ time windows 
overlap is illustrated in fig. [7[ 

A.2 Time-domain method 

After having downsampled^ the data, the peaks of the sinusoidal temperature streams at the two 
points are found numerically. From each consecutive pair of peaks an estimate of frequency v, 
damping y and phase shift Ay? can be easily derived. Assuming that we have A'^ pairs, the best 
values for v, y and Aip are simply the average of the A'^ estimates. Basic eiTor propagation is used 
to derive estimates for the eiTors associated with each quantity. 

This method has the advantage of providing a straightforward way to predict the order of mag- 
nitude of the expected error on the measurement of y and before actually doing the calculations. 
As an example, we can derive the error on the y parameter considering the accuracy of the focal 
plane thermometers: if this accuracy is 20 mK, then the fluctuation amplitude (difference between 
two temperatures) has an error of 

V2x(20mK)2 _ 20mK 
^/2N y/N ' 

where A'^ is the number of periods and is typically between 15 and 30. Since the absolute level of 
the temperature during the RAA tests was between 20 and 30 K, this means that the relative error 
on the amplitude is of the order of 10"^: this is roughly the same as the expected error for y, since 
this quantity is a ratio between two amplitudes. These are indeed the eiTors we get on this quantity 
from the analysis of the Planck/LFI RAA tests (see table 

^We chose not to divide this error by VA' because the A' estimates are not truly independent, as they are calculated 
over time windows which partially overlap. 

^The frequency of the downsampled data is chosen iteratively, with the goal of optimizing the relative error on the 
estimates for y and Aip. 




Figure 8. Temperature profile for a two-body system in equilibrium extending from x-Otox-L and with 
a contact resistance Rc at x = 1 (see eq. |B.5| and |B.6| ). 



A.3 Fitting method 

Temperature profiles are fitted numerically with a function of the form T sm{2nvt + 9) + p{t), where 
p{t) is a polynomial which models slow temperature drifts. From this fit the values for y and A<^ as 
well as their errors can be easily estimated. 



B. Analytical Models of Contact Resistances 



In this appendix we discuss two simple analytical models that show the effects due to a non-ideal 
thermal contact between two bodies (hereafter named "body 1" and "body 2"). We only use di- 
mensionless quantities for the sake of simplicity. The bodies are supposed to be monodimensional 
along the x direction and connected at x - 1, and their thermal properties (heat capacity, density 
and thermal conductivity) do not depend on x nor on time t. The governing equations for the 
temperatures Ti and T2 in the two bodies are: 

dlji{x,t) = d,Ti{x,t), (B.l) 

ddlj2{x,t) = dtT2{x,t), (B.2) 

-Rcd.Tiiht) = Ti{\,t)-T2(\,t), (B.3) 

-kRcd,T2{\,t) = Ti{l,t)-T2{l,t), (B.4) 

where d and k are the thermal diffusivity and conductivity of the second body with respect to the 
first and Rc is the contact resistance. Equation |B.1| and |B.2| are the heat equation expressed for each 
of the two bodies, equations |B.3| and |B.4| use the definition of contact resistance, which is the ratio 
per unit surface of the temperature jump at the boundary over the heat flux q = -kdxT . 



B.l Static case 

A very simple solution for the set of equations |B.1HB.4| is found in the steady-state case {dfTi - 
5(7^2 - 0). If body 1 extends from x-Qtox- 1, and body 2 has a length L, then from the boundary 
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conditions Ti\^^q - and T2\x=\+l - 1 we derive the following solutions: 

k{l +Rc) + L 

Tiix) = 1 - (B.6) 
k(l +Rc) + L 

where Ti{x) is defined for x e [0, 1) and T2{x) for x e (1, 1 + L], i.e. there is a sudden jump in the 
temperature across the boundary, and the temperature varies linearly in the two bodies with a slope 
which depends on the thermal conductivity of the body. This solution is plotted in fig. ^ 

B.2 Dynamic case 

In analogy with the Planck/LFI tests discussed in this paper, we consider here a sinusoidal temper- 
ature fluctuation propagating through the two bodies. The fluctuation is applied to body 1 at .s; = 0, 
and we assume that body 2 is infinitely long^: 

Ti{x^O,t) = Tae~^''"', (B.7) 
lim \T2{x,t)\ = 0, (B.8) 

.V— >+oo 

with Ta the temperature fluctuation amplitude. 
We write the generic solution Ti /2 as 

T„2(x,t)^Tagi/2{x)e-^''''". 

Under this hypothesis, |ragi/2(.x)| is the fluctuation amplitude of the temperature at point x and 
arg^i/2(ji:) is the phase of the fluctuation. Equations from |B.1| to |B.4| rewrite as follows: 

g'{{x) = 2mvg,{x), (B.9) 

dg'{{x) = 2mvg2{x), (B.IO) 

-Rcg[W = giW-giW, (B.ll) 

-kRcg'2{l) = gi{l)-g2W. (B.12) 

Putting a = ^Jkv, ^ -Rc ^^firv, we get the solutions 

glix) - ciexp(a{l+i)x^ + {l-ci)exp(^-a{l + i)x^, (B.13) 

glix) = C2exp|--^(1+;>J, (B.14) 



where 



{i + i)^k-k+Vd 

Cl = 7 (B.15) 

(1 + i)^k -k + e2«(i+;) ((1 + i)^k + k+yfd)+y/d 

2V5exp(a(l.,)(l.^)) 

C2 = } (B.16) 

(1 + i)^k -k + f-Mi+O ((1 + i)^k + k+yfd)+yjd 



^The boundary condition for Ti has a minus sign in tiie exponential in analogy with the standard way to write plane 
waves, expi(kx- (jjt). 
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Figure 9. Fluctuation amplitude (Igi/iC-^)!) and phase (arggi/2(x)) for a sinusoidal fluctuation propagating 
from X = and encountering a discontinuity with contact resistance Rc at jc = 1 (equations |B.13| and |B.14| ). 
Three values of Rc have been used to produce the plots. The values used for the parameters are a = \,k-2, 
d=l/2. 



The amplitude and phase of the fluctuation in the two bodies as a function of x are shown in fig. | 
At the boundary (x = 1) the ratio between giix) and g2(x) is the complex number 



lim 



glix) 



1 



lim^^i- glix) y/d/{k^) + \ + f 

whose absolute value 7rc and argument Ai^rc are 



(B.17) 



1 



\ {^/{ko+lf+l 

Av^Rc = -arctan 



I yfd/ik^+l 



(B.18) 
(B.19) 



The value of 7rc is the "amplitude drop" at x = I (left side of fig. 0) and the value of ^rc is the 
phase shift at the same point (right side of the same figure). Therefore, fig. shows that a non-ideal 
thermal contact between two bodies has the eff'ect of damping the amplitude of thermal fluctuations 
and inducing a phase delay in the propagating wave. Both eff'ects increase for large values of 
(i.e. if Rc +00 then — > and Aip — > -n/4). 

It is easy from these formulae to obtain the expressions with the proper dimensions. If the two 
bodies have thermal difi'usivity coefficients Di and D2 and thermal conductivities k\ and k2, and 
the discontinuity is at .x = L, it is enough to use the following substitutions in equations from p.l3| 
to EM. 

^ ^2 , k2 [nv . „ . [ttv 

d= — , k= — , a=. — , t=RcA. — , 
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where A is the area of the boundary surface at x = L. Moreover, in equations |B.15| and |B.16| the 
substitution 




(B.20) 



must be made in the exponentials. 
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Figure 10. Comparison between the estimates of the numerical model and the experimental results for the 
y function (damping). X axis: frequency in mHz, Y axis: value of y (pure number). Error bars have been 
magnified by a factor 50 in order to make them visible. 
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Figure 11. Comparison between the estimates of the numerical model and the experimental results for the 
Atp function (phase shift). X axis: frequency in mHz, Y axis: value of A(p (radians). Error bars have been 
magnified by a factor 10 in order to make them visible. 
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Table 4. Estimated value of v, y and A(p for each of the three frequencies tested in the laboratory. This 
table only reports the best estimate among the ones obtained with the three methods (time-domain, Fourier, 
fitting) described in the text. The detailed results for y and Aip are reported in tables § and 

Sensor v [mHz] y [rad] 

TSIL 0.1995±0.0008 0.4168±0.0008 0.551±0.011 

0.5997±0.0009 0.2756±0.0004 0.7 11 ±0.008 

1.3999±0.0008 0.2554±0.0003 0.960±0.019 

TS3L 0.1995±0.0008 0.5023±0.0010 0.584+0.012 

0.5995±0.0012 0.3125±0.0005 0.729±0.008 

1.4002±0.0008 0.2231±0.0003 0.880±0.010 

TS4L 0.1995±0.0008 0.4853±0.0010 0.632±0.013 

0.5996±0.0010 0.2919±0.0004 0.844±0.009 

1.4001±0.0011 0.1931±0.0002 1.100±0.010 

TS6L 0.1995±0.0008 0.5138±0.0010 0.550±0.011 

0.5996±0.0010 0.3452±0.0005 0.708±0.008 

1.4001 ±0.0004 0.2576±0.0003 0.970±0.010 

TSIR 0.1995±0.0008 0.5463±0.0011 0.733±0.013 

0.5995±0.0013 0.2971±0.0006 1.059±0.010 

1.3991±0.0025 0.1674±0.0002 1.371±0.022 

TS2R 0.1995±0.0008 0.5658±0.0011 0.839±0.016 

0.5994±0.0015 0.3124±0.0007 1.292±0.012 

1.3996±0.0017 0.1583±0.0002 1.841±0.022 

TS3R 0.1995±0.0008 0.5795±0.0012 0.860±0.016 

0.5991±0.0021 0.2918±0.0009 1.385±0.011 

1.4002±0.0016 0.1317±0.0002 1.990±0.010 

TS4R 0.1995±0.0008 0.4725±0.0009 0.677±0.013 

0.5996±0.0011 0.2666±0.0005 0.933±0.009 

1.3998±0.0011 0.1686±0.0002 1.160±0.022 

TS5R 0.1995±0.0008 0.4746±0.0010 0.668±0.013 

0.5997±0.0009 0.3411 ±0.0005 0.934±0.009 

1.3999±0.0008 0.2194±0.0003 1.276±0.022 

TS6R 0.1994±0.0008 0.4954±0.0009 0.918±0.017 

0.5994±0.0014 0.2433±0.0005 1.510±0.012 

1.3990±0.0026 0.1202±0.0002 2.277±0.018 
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Table 5. Estimated value of y for each of the three frequencies tested in the laboratory. AH values are in 
radians. The results provided by each analysis method are reported here. 



Sensor 



TSIL 



TS2L 



TS3L 



TS4L 



TS5L 



TS6L 



TSIR 



TS2R 



TS3R 



TS4R 



TS5R 



TS6R 



Fourier 



Time-domain 



Fit 



0.417±0.002 
0.275±0.002 
0.256±0.004 



0.416±0.004 
0.273±0.003 
0.257±0.003 



0.4168±0.0008 
0.2756±0.0004 
0.2554±0.0003 



0.527±0.002 

0.432±0.001 
0.478±0.004 



0.526±0.004 

0.433±0.003 
0.477±0.004 



0.5282±0.0008 

0.4346±0.0005 
0.4790±0.0005 



0.502±0.003 
0.3 11 ±0.002 
0.222±0.004 



0.500±0.005 
0.312±0.004 
0.225±0.004 



0.5023±0.0010 
0.3125±0.0005 
0.2231 ±0.0003 



0.486±0.003 

0.291 ±0.002 
0.193±0.004 



0.483±0.005 

0.293±0.003 
0.192±0.004 



0.4853±0.0010 

0.2919±0.0004 
0.1931 ±0.0002 



0.536±0.002 
0.448 ±0.001 
0.532±0.004 



0.538±0.003 
0.451 ±0.002 
0.528±0.003 



0.5378±0.0008 
0.4479±0.0005 
0.5301±0.0006 



0.515±0.003 

0.344+0.002 
0.259+0.004 



0.507±0.004 

0.345+0.003 
0.262+0.002 



0.5138±0.0010 

0.3452+0.0005 
0.2576±0.0003 



0.547±0.003 
0.298±0.003 
0.168±0.004 



0.545±0.005 
0.294±0.004 
0.171±0.003 



0.5463±0.0011 
0.2971 ±0.0006 
0.1674±0.0002 



0.569±0.004 
0.313±0.003 
0.158±0.003 



0.564±0.006 
0.313±0.005 
0.159±0.004 



0.5658±0.0011 
0.3124±0.0007 
0.1583±0.0002 



0.582±0.004 
0.293±0.003 
0.130±0.003 



0.576±0.006 
0.291 ±0.004 
0.131±0.004 



0.5795±0.0012 
0.2918±0.0009 
0.1317±0.0002 



0.474±0.003 
0.264±0.002 
0.169±0.003 



0.471 ±0.003 
0.265±0.002 

0.170±0.002 



0.4725±0.0009 
0.2666±0.0005 

0.1686±0.0002 



0.476±0.003 
0.340±0.003 
0.219±0.004 



0.473±0.003 
0.340±0.003 
0.222±0.002 



0.4746±0.0010 
0.3411 ±0.0005 
0.2194±0.0003 



0.498±0.003 
0.246±0.002 
0.119±0.002 



0.493±0.003 
0.243±0.003 
0.118±0.003 



0.4954±0.0009 
0.2433±0.0005 
0.1202±0.0002 
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Table 6. Estimated value of for each of the three frequencies tested in the laboratory. All values are in 
radians. The results provided by each analysis method are reported here. Note how the fitting method shows 
the largest discrepancies, although its error bars are always compatible with the other methods. 



Sensor 


Fourier 


Time-domain 


Fit 


TSIL 


0.551+0.011 


0.548+0.042 


0.73+0.18 




0.7 11 ±0.008 


0.715+0.028 


0.61+0.20 




0.960±0.019 


0.974+0.044 


0.94+0.03 


TS2L 


0.315±0.006 


0.306+0.034 


0.47+0.15 




0.311+0.004 


0.309+0.023 


0.28+0.12 




0.438+0.01 1 


0.440+0.034 


0.43+0.01 


TS3L 


0.584±0.012 


0.577+0.044 


0.77+0.19 




0.729+0.008 


0.732+0.032 


0.57+0.25 




0.885+0.018 


0.879+0.048 


0.88+0.01 


TS4L 


0.632+0.013 


0.617+0.045 


0.83+0.20 




0.844+0.009 


0.856+0.034 


0.73+0.21 




1.129+0.021 


1.114+0.052 


1.10+0.01 


TS5L 


0.301+0.006 


0.298+0.026 


0.46+0.16 




0.290+0.004 


0.295+0.020 


0.27+0.12 




0.428+0.010 


0.432+0.032 


0.42+0.02 


TS6L 


0.550+0.011 


0.544+0.032 


0.73+0.18 




0.708+0.008 


0.708+0.024 


0.60+0.21 




0.983+0.019 


0.990+0.042 


0.97+0.01 


TSIR 


0.733+0.013 


0.727+0.039 


0.93+0.20 




1.059+0.010 


1.059+0.033 


0.89+0.26 




1.371+0.022 


1.376+0.059 


1.27+0.09 


TS2R 


0.839+0.016 


0.828+0.045 


1.05+0.21 




1.292+0.012 


1.301+0.036 


1.09+0.30 




1.841+0.022 


1.802+0.067 


1.76+0.05 


TS3R 


0.860+0.016 


0.849+0.043 


1.08+0.22 




1.385+0.011 


1.389+0.035 


1.03+0.44 




1.982+0.024 


1.988+0.075 


1.99+0.01 


TS4R 


0.677+0.013 


0.660+0.031 


0.87+0.19 




0.933+0.009 


0.930+0.026 


0.79+0.24 




1.160+0.022 


1.178+0.044 


1.13+0.03 


TS5R 


0.668+0.013 


0.660+0.030 


0.87+0.20 




0.934+0.009 


0.936+0.025 


0.85+0.18 




1.276+0.022 


1.260+0.042 


1.25+0.03 


TS6R 


0.918+0.017 


0.903+0.033 


1.14+0.21 




1.510+0.012 


1.523+0.029 


1.32+0.29 




2.277+0.018 


2.271+0.060 


2.18+0.10 
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